Comments on "Modeling Galaxy Halos Using Dark Matter with Pressure" 



o 
o 

(N 

>>: 



Kung-Yi Sui'=0 and Pisin Cheni'^.s.'jt] 
1. Department of Physics, National Taiwan University, Taipei, Taiwan 10617 
2. Graduate Institute of Astrophysics, National Taiwan University, Taipei, Taiwan 10617 
3. Leung Center for Cosmology and Particle Astrophysics, 

National Taiwan University, Taipei, Taiwan 10617 
4- Kavli Institute for Particle Astrophysics and Cosmology, 
SLAC National Accelerator Laboratory, Stanford University, Stanford, CA 94305, U.S.A. 

(Dated: May 14, 2009) 

We comment on the calculational mistake in the paper "Modeling galaxy halos using dark matter 
with pressure" by Somnath Bharadwaj and Sayan Kar. The authors made a mistake while calculat- 
ing the metric, which led to an overestimate of the deflection angle of light passing through the halos 
for — 1 < TOr < —0.5 and an underestimate of the deflection angle for —0.5 < Wr < 0. In addition, 
solution for w,. > should not exist. Although the Bharadwaj-Kar solution should be corrected, it 
appears that the characteristics of the deflection angle under the supposed non-conventional non- 
ideal fluid equation of state for the dark matter halo remain sensitive to the impact parameter and 
may be veriflable through observations. 

PACS numbers; 98.35.Gi, 95.30.Sf 
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In the paper "Modeling galaxy halos using dark matter 
with pressure" [ij , Somnath Bharadwaj and Sayan Kar 
inspected the possibility that the dark matter could be a 
non-ideal fluid with a significant pressure. They showed 
that such a non-conventional notion of dark matter could 
be verified through the observed galactic rotation curves 
and the modification on the gravitational lensing deflec- 
tion angle. This is interesting and different from the con- 
ventional approach to dark matter. Such a dark matter 
cannot be accommodated by the Newtonian gravity and 
it has to resort to general relativity. They model the 
galactic dark matter halo based on the following assump- 
tions: 

1. The metric inside the halo is spherically symmetric 
and the proper time interval dr is 



(1) 



and it reduces to Schwarzschild metric outside the 
boundary, r = i?, of the halo. 

2. The gravitation is weak. 

3. The pressure inside the halo is anisotropic, which 
means that the radial pressure Pr is different from 
the tangential pressure Pt. 

4. The equation of state is either P^ = Wrpc^ or Pt = 
Wtpc^, where Wr and Wt are constants and p is the 
energy density. 



A. Metric 

The task is to determine goo and grr, or <I>(r) and A(r), 
in the metric. and therefore goo, can be fully de- 

termined by the flat rotation curve. A(r), on the other 
hand, can be determined by the equation of state with 
the aid of the Einstein equations. After substituting the 
$(r), the Einstein equations lead to three equations for 
A(r): 
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where Vc ~ 200 km/s is the constant rotation velocity 
of the halo associated with the flat rotation curve of our 
galaxy. Using the above metric, one flnds that the de- 
flection angle depends strongly on the equation of state. 
Accordingly with the aid of the observed deflection angle, 
the equation of state can be determined. Conversely, if 
the equation of state of the dark matter halo is known, 
one can in principle test the general relativity at galactic 
scale. 

While we agree with the framework of their approach, 
we find some mistakes in their calculation of the metric 
under the equation of state Pr = Wrpc^, which result in 
an overestimate of the defiection angle of light passing 
through the halo for —1 < Wr < —0.5 and an underes- 
timate of it for —0.5 < Wr < 0. In addition, solution 
for Wr > should not exist. On the other hand, we 
found no mistake in their calculation under the condi- 
tion Pt = Wtpc^, and the corresponding conclusions for 
that case remain unchanged. 

Their original solutions and our corrected expressions 
of the metric, the density and the pressure are displayed 
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TABLE I. Bharadwaj-Kar solution (r < R) 



Wr 


A(r) = AjvX p{r)=pNX 


Pt{r) = {c^/2)pN{r)x 


— l<Wr<O0 
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1 - ln{r/R) ln{r/R) 


1 


TABLE II. The corrected solution (r < R) 


Wr 


A(r) = XnX p{r) = pnX 


A(r) = (cV2)p^(r)x 


-KWr <Q 






-1 


1 - \n{r/R) ln(r/R) 


1 



in Table 1 and Tabic 2, respectively, for comparison. 
Note that Ajv = {vc/cf and pjv = v'^/AnGr'^. 

It can be seen that the main difference between the two 
results lies in the second term of A, where the exponent 
should be —(1 + Wr)/wr instead of — + Wr)- As a 
consequence, the acceptable range of Wr should be [— 1 , 0] 
rather than [— 1 , oo) , so as to ensure the positive definite- 
ness of the energy density. It should be evident that the 
two solutions are identical when Wr = —0.5. Even so, 
several qualitative features of the Bharadwaj-Kar solu- 



tion still hold in our results. For example, for Pj. 

-9oo 



the metric elements satisfy the relation gn = —Pnn inside 
and outside the halo in both calculations. 

Bharadwaj and Kar commented in their paper that 
Pr = —WrP(? is the only condition for p = at r = 
R. They also mentioned that in general the pressure 
and the energy cannot be matched at the halo boundary 
and would therefore be discontinuous. However we see 
from the corrected solution that the situation is not so 
bleak. For the equation of state Pr = Wrpc^, we find 
that both the radial pressure and the energy density are 
and continuous at the boundary, although the tangential 
pressure remains discontinuous. As for the equation of 
state Pt = wtpc? , the radial pressure is and continuous 
at the boundary but the energy density and tangential 
pressure are not. 
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FIG. 1: The deflection angle under Pr = Wrpc^ based on the 
Bharadwaj-Kar solution. 
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B. Deflection Angle 



FIG. 2: The corrected deflection angle under Pr = Wrpc^ ■ 



The general formula for the deflection angle in the weak 
field limit, i.e., Eq(15) in the original paper, is 



(3) 



where 5 is the deflection angle evaluated in units of 
{vc/cf, y = u + a{y), a{y) = [($ + X)/b^ - Xy^]/y and 
u = l/r. b is the impact parameter that equals to L/E, 
where L = r^{d^/dT) and E = e^"^ c{dt / dr) are the con- 
stants of motion of the system. For a given equation 
of state, the deflection angle can be obtained based on 
the above formula. Fig.l and Fig. 2 display the author's 



original and our corrected results, respectively, of the de- 
flection angle as a function of b/R under the equation of 
state Pr = WrPC^- 

For the case Pr = Wrpc^, several features of the deflec- 
tion angle remain unchanged in both solutions. When 
Wr < 0, the deflection angle is always larger than that 
in the case where there is no pressure. Furthermore, the 
deflection angle increases as b/R decreases, and it ap- 
proaches a constant that is larger than Sn = 27r, which is 
the deflection angle for an infinite halo with no pressure. 
The main difference, however, is that the corrected de- 
flection angle converges more slowly to a smaller value for 
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Wr < —0.5, while faster to a larger value for Wr > —0.5. 
We note that in Fig. 2 the deflection angle for Wr — —0.9 
in our calculation does converge, though very slowly. In 
addition, we see from the figures that there exists a "turn- 
ing point " for some of the curves. In the Bharadwaj-Kar 
solution, the turning point is more visible for larger val- 
ues of Wr, while it is more visible for smaller values in 
our calculation. 

In conclusion, our corrected solution for the case Pr — 



Wrpc^ leads to smoother energy and pressure profiles at 
the halo boundary. In addition, the corresponding deflec- 
tion angle converges more slowly to a smaller value for 
— 1 < Wr < —0.5, while it converges more quickly to a 
larger value for —0.5 < Wr < 0. The good news is, under 
the new solution the deflection angle remains sensitive to 
the equation of state. So it may still be possible to verify 
the proposed non-conventional, non-ideal fluid equation 
of state for the dark matter halo through observations 0. 
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